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摘 要 : 干旱 深刻 影响 着 巴基斯坦 的 水 资源 、 粮 食 产 量 和 社会 经 济 格局 ,分 析 当 地 的 干旱 特征 及 其 
风险 对 其 社会 发 展 具有 重要 意义 ,但 当前 鲜 有 学 者 从 长 期 视角 分 析 整 个 巴基斯坦 的 干旱 特征 或 评 
估 该 地 区 的 干旱 灾害 风险 。 通 过 采用 1982 一 2013 年 第 三 代 归 一 化 植被 指数 (Global inventory mod- 
elling and mapping studies normalized difference vegetation index ,GIMMS NDVI3g) 数 据 集 构建 植被 状 
态 指 数 (Vegetation condition index, VCI) ,分 析 了 32a 书 基 斯 坦 干旱 范围 和 频率 ;采用 由 干旱 频率 、 
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风险 进行 了 评估 。 结 果 表 明 :(1) 研究 时 段 内 ,各 月 多 年 平均 干旱 范围 存在 差异 ,1 月 和 9 一 12 HT 


旱 范 围 较 大 (面积 比例 大 于 40% );5 一 8 月 干旱 范围 较 小 (天 


水 减少 影响 ,干旱 
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积 比例 小 于 30%);5 一 6 月 受 升温 和 降 


RA LF HH (2) 研究 区 5 一 8 月 干旱 频率 较 低 ,西部 山区 林地 、 草 地 9 月 干旱 
原 1 月 和 9 一 12 月 干旱 频率 较 高 ,。(3) 巴基斯坦 干旱 风险 空间 格局 主要 由 干旱 频 


率 、 作 物产 量 、 大 牲畜 比例 、 土 壤 持 水 能 力 决定 ,其 中 干旱 频率 对 区 域 干旱 风险 的 影响 最 大 。 在 牧 
畜 比 例 和 土壤 持 水 性 的 分 别 作用 下 ,巴基斯坦 北部 山区 向 平原 过 渡 的 地 带 和 东南 部 自然 植被 覆盖 
地 区 干旱 风险 较 高 ; 受 灌溉 影响 ,地 形 平坦 的 河谷 地 区 干旱 风 险 较 低 ;北部 山区 自然 植被 覆盖 度 


高 ,于 旱 风险 也 较 低 。 研 究 结果 有 益 于 巴基斯坦 灾害 风险 管理 。 
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干旱 作为 众多 自然 灾害 的 1 种 ,与 火山 地震、 
洪涝 飓风 等 典型 突 发 性 自然 灾害 相 比 ,其 发 生 具 
有 渐进 性 和 潜在 性 , 且 持 续 时 间 长 "。 严 重 且 持 续 
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的 干旱 将 造成 士 地 荒漠 化 .河水 断 流 等 现象 ,从 而 
导致 农业 减产 ,制约 社会 经 济 发 展 ”1。 近 年 来 , 巴 
基 斯 坦 气 温 上 升 ,季风 降水 量 减少 ,作物 需 水 不 断 
增长 ,作物 需 水 量 已 经 超过 了 有 效 降 水 量 ,导致 干 
旱灾 害 频 发 ,严重 制约 当地 农业 的 发 展 。 在 全 球 
气候 变 暖 背景 下 ,对 巴基斯坦 干旱 特征 及 其 风险 进 
行 分 析 , 具 有 重要 的 现实 意义 ,将 有 益 于 该 地 区 灾 
害 风险 管理 .农业 生产 、 社 会 经 济 可 持续 发 展 和 生 
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不 同 地 区 的 干旱 展开 了 人 研究 。 如 Lu 等 ”基于 
1981—2013 年 植被 覆盖 指数 (Clobal inventory mod- 
elling and mapping studies advanced very high resolu- 
tion radiometer normalized difference vegetation index , 
GIMMS AVHRR NDVI) 长 时 间 序 列 数据 集 ,综合 北 
美 区 域 再 分 析 (North American regional reanalysis , 
NARR ) 的 地 表 温 度 、 降 水 .土壤 含水 量 等 数据 ,开发 
了 综合 干旱 指数 (Integrated scaled drought index , IS- 
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DD ,对 美国 的 干旱 分 布 情况 进行 了 分 析 ;Zhang 等 * 
利用 印度 1981—2013 年 月 尺度 GIMMS NDVI, 获得 
了 植被 状态 指数 (VCI) 时 间 序 列 ,研究 了 印度 地 区 
干旱 持续 时 间 频率 .范围 和 严重 程度 的 发 展 趋势 ， 
并 建立 了 作物 产量 异常 与 干旱 之 间 的 关系 。 也 有 学 
者 将 遥感 运用 到 巴基斯坦 干旱 研究 中 ,如 Shah 等 ” 
利用 1998 一 2014 年 的 MODIS NDVI 数 据 生成 干旱 
指数 (Standardized vegetation index, SVI) ,对 巴 基 斯 
JH Tharparkar 地 区 的 干旱 严重 程度 进行 了 时 空 4 
析 。 除 遥感 数据 外 ,Adnan 等 "还 基于 全 球 气 候 中 
心 的 降水 数据 以 及 气候 预测 中 心 1951—2010 年间 
的 土壤 湿度 数据 ,利用 区 域 干旱 识别 模型 分 别 对 巴 
基 斯 坦 Sindh 省 3 个 月 .6 个 月 .9 个 月 .12 个 月 .24 个 
月 的 标准 化 降水 指数 (Standardized precipitation in- 
dex ,SPI) 进 行 了 计算 ,利用 SPI 分 析 了 该 地 区 的 干旱 
特征 。 巴 基 斯 坦 纬度 跨度 大 ,气候 地 形 条 件 不 一 ， 
区 域 空间 差异 大 ,但 当前 大 部 分 工作 侧重 于 巴 共 斯 
坦 部 分 地 区 的 干旱 监测 以 及 干旱 时 空 分 析 , 缺 乏 对 
整个 巴基斯坦 干旱 的 长 期 视角 研究 分 析 , 且 没有 从 
灾害 风险 角度 来 评估 巴基斯坦 干旱 灾害 风险 。 

鉴于 此 ,本 文采 用 1982 一 2013 年 GIMMS ND- 
VIl3g 长 时 间 序 列 数据 集 获取 VCI, 在 利用 VCI 对 巴 
基 斯 坦 干 旱 特征 进行 分 析 的 基础 上 ,通过 干旱 频 
率 .作物 种 植 面 积 比 例 灌溉 比例 作物 产量 .田间 
持 水 量 、 人 口 密度 、 大 冀 牲 比例 等 资料 ,构建 巴 基 斯 
JH T fe E .脆弱 性 .暴露 度 以 及 人 类 防 灾 减 灾 
能 力 等 干旱 风险 性 指标 ,建立 干旱 风险 评估 模型 ， 
对 巴基斯坦 全 域 进行 干旱 风险 评估 ,以 期 为 巴 基 斯 
坦 自然 灾害 风险 防 控 与 减灾 对 策 的 制定 提供 科学 
依据 。 


1 研究 区 概况 


巴基斯坦 地 处 南亚 ,地 理 范围 为 60°~80°E、23°~ 
37°N, 总 面积 88x10' km? (包括 巴 控 克什米尔 地 
区 )。 从 来 源 于 美国 地 质 测量 局 (USGS) 的 GTO- 
P030 数字 高 程 图 (图 1a) 可 以 看 出 ,巴基斯坦 多 山 
地 丘陵 ,北部 为 山地 、 高 原 ,地 势 较 高 ,西部 为 山 
地 ,中 部 、 南 部 为 平原 。 由 于 地 貌 条 件 复杂 ,南北 纬 
度 跨度 大 ,巴基斯坦 气候 种 类 较 多 ,北部 为 亚热带 
气候 和 高 山 气 候 , 西 南部 为 热带 沙漠 气候 ,东南 部 
为 热带 季风 气候 。 巴 基 斯 坦 区 域 平均 年 降水 量 不 
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足 250 mm,1/4 地 区 降水 量 在 120 mm 以 下 ;特别 是 
夏季 , 受 副 高 控制 ,气候 炎热 干燥 ,降水 较 少 。 根 据 
ESA CCI 2010 土 地 覆盖 数据 (http://maps.elie.ucl.ac. 
be/CCI/viewer/ ) ,巴基斯坦 土地 覆盖 类 型 (图 1b) 主 
要 包括 灌溉 农田 草地、 林地 荒漠 等 ,面积 比例 分 
别 为 22% 24% 12% , 27796 ,农田 主要 位 于 东部 、 南 
部 河谷 地 区 ,荒漠 主要 位 于 西部 Balochistan 省 境 
内 。 巴 基 斯 坦 土壤 质地 类 型 …( 图 le) 主要 包括 砂 
土壤 砂 土 , 粉 壤 壤土 ` 黏 壤 、 粉 黏 壤 等 ,农田 地 区 
以 壤土 BD RE RRN E RN HE e UHR UA 
砂 士 和 壤 砂 士 为 主 。 巴 基 斯 坦 人 口 自 然 增长 率 为 
1.55% ,根据 文献 图 集中 1998 年 人 口 分 布 和 该 增 
长 率 , 获 得 2010 年 巴基斯坦 人 口 分 布 (图 1d)。 从 图 
中 可 以 看 出 ,巴基斯坦 人 口 主要 分 布 于 中 东部 的 
Sindh 省、Punjab 省 、Khyber-Pakhtunkhwa 省 等 地 。 


2 数据 与 方法 


本 研究 利用 由 遥感 植被 指数 得 到 的 干旱 指数 
分 析 巴 基 斯 坦 干旱 ,并 在 评估 再 分 析 气 象 数 据 的 基 
础 上 ,利用 再 分 析 数 据 研究 干旱 驱动 因子 。 此 外 ， 
本 研究 将 利用 田间 持 水 量 灌溉 面积 比例 人口 密 
BE .大 牲畜 比例 、 作 物种 植 面 积 比例 .干旱 频率 、 作 
物 (小 麦 .玉米 水稻 .棉花 ) 产 量 等 数据 进行 干旱 风 
险 性 指标 因子 计算 ,评估 巴基斯坦 干旱 风险 。 

2.1 数据 

2.1.1 遥感 植被 指数 NDVI 数 据 能 够 很 好 的 反映 
植被 的 季节 和 年 变化 ” ,在 全 球 和 区 域 植被 变化 及 
干旱 监测 中 得 到 广泛 应 用 “”""。 本 人 研究 利用 1982 一 
2013 年 GIMMS NDVI3g 长 时 间 序 列 数据 集 分 析 巴 
基 斯 坦 干旱 ,该 数据 时 间 分 辨 率 为 15 q ,空间 分 辨 率 
为 0.0833°(http://ecocast.arc.nasa.gov/data/pub/gimms/ 
3g.v0/) ,采用 Lambert 等 积 方位 投影 ,利用 最 近邻 法 
内 捅 重 采 样 至 8 km。 再 利用 最 大 值 合成 法 得 到 月 
NDVL, 进 一 步 降低 大 气 残 余 效 应 及 双向 传输 效应 的 
影响 *。 

2.1.2 气象 数据 本 人 研究 中 的 气象 数据 包括 再 分 析 
资料 与 气象 站 实测 资料 。 再 分 析 资 料 来 源 于 由 欧 
洲 中 期 数值 预报 中 心 再 分 析 数 据 (ERA-Interim ) 经 
海拔 和 月 尺度 实测 数据 矫正 ”得 到 的 时 间 分 辩 率 
为 1 d, 空 间 分 辨 率 为 0." 的 全 球 实测 气候 数据 集 


(http://www.waterandclimatechange.eu/about/watch- forc- 
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Fig. 1 General situation of Pakistan 


ing-data-20th-century ) ,实测 资料 (来 源 于 NOAA 观 
测 台 站 数据 ,http://ftp.ncdcnoaa.gov/pub/data/gsod/) 
包括 研究 时 段 内 资料 相对 完整 的 巴基斯坦 10 个 气 
象 站 点 (图 2a 和 2c) 的 气温 和 降水 数据 。 在 利用 
ERA-Interim 再 分 析 数 据 研究 巴基斯坦 干旱 之 前 ,本 
研究 将 采用 站 点 周围 栅 格 面积 加 权 平 均 法 评估 
数据 集 在 该 地 区 的 精度 。 利 用 实测 资料 对 再 分 析 
数据 的 评估 表明 ,再 分 析 气 温和 降水 能 够 反映 各 站 
点 实际 月 均 温 (图 2a) 和 月 降水 (图 2c) 的 变化 (相关 
系数 大 于 0.48 ,通过 了 5% 水 平 显著 性 检验 ) , 且 对 气 
温 的 重 现 性 能 强 于 降水 ,再 分 析 数 据 也 能 够 反映 区 
域 平均 月 均 温 (图 2b) 和 月 降水 (图 2d) 的 变化 (决定 
系数 大 于 0.9)。 

2.1.3 农业 数据 为 了 分 析 人 研究 区 内 的 干旱 风险 ， 
还 利用 了 作物 种 植 面积 比例 \ 作 物 综合 产量 (作物 
种 植 面积 比例 加 权 平 均 作 物 单 产 ) 洪波 面积 比例 、 
大 冀 牲 比例 等 数据 (来 源 于 巴基斯坦 国家 统计 局 及 


文献 图 集 "")。 从 图 3 可 以 看 出 ,巴基斯坦 作物 种 植 
区 主要 位 于 东部 的 Punjab 省 、Sindh 省 以 及 首都 Is- 
lamabad 附近。 巴基斯坦 主要 种 植 棉花 、 小麦 .水 
稻 、 玉 米 等 ,棉花 -小 麦 为 主要 作物 种 植 模式 。 作 物 
产量 较 高 地 区 主要 位 于 东南 部 的 Punjab 省 和 Sindh 
省 ,两 省 作物 综合 产量 分 别 为 2890 kg'hm” 和 3889 
kg.hm ,其 中 Punjab 省 小 麦 产量 为 2694 kg:hm ,7K 
稻 产 量 为 1842 kg.hm ,玉米 产量 为 4916 kg* hm”, 
棉花 产量 为 3935 kg- hm? Sindh 省 小 麦 水稻、 玉米 、 
棉花 作物 产量 分 别 为 3432 kg -hm 7,3458 kg hm’, 
629 kg.hm? 5304 kg.hm?。 巴 基 斯 坦 灌 溉 面积 比例 
较 大 (>80% ) 的 地 区 主要 位 于 Punjab 中 东部 。 由 于 
西部 地 区 为 草地 和 荒漠 ,适合 小 牲畜 ,东部 地 区 为 
林地 和 农田 ,适合 大 牲畜 ( 牛 ), 因 此 小 牲畜 ( 羊 ) 主 
要 位 于 西部 地 区 ,大 畜生 主要 位 于 东部 及 北部 地 
区 。 这 些 数 据 主要 用 于 干旱 胸 弱 性 .干旱 暴露 度 、 
人 类 防 灾 减灾 能 力 等 指标 的 计算 ,从 而 构建 巴 基 斯 
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图 2 再 分 析 月 均 温和 降水 评估 


Fig.2 Evaluation of ERA-Interim temperature and precipitation at the monthly time scale 
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图 3 巴基斯坦 农业 示意 图 


Fig. 3 Pakistan agricultural information 
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坦 的 干旱 灾害 风险 性 指数 。 
2.2 方法 
2.2.1 和 干旱 监测 指数 ”在 众多 干旱 监测 指数 中 , 植 
被 状态 指数 (VCI) 能 够 反映 植被 生长 状况 ,在 干旱 
监测 中 应 用 较为 广泛 2 ,因此 本 研究 将 利用 VCI 
对 巴基斯坦 干旱 进行 监测 。VCI 计 算 公 式 如 下 
所 示 : 

yora NOV (D 

i- NDVI,,, - NDVI,,, 

式 中 :NDVJ 为 某 年 时 期 的 NDVI 值 ;NDVI、NDVI 
分 别 为 1982—2013 年 时 期 7 的 NDVI 最 大 值 和 最 小 
值 。 

VCI 指 数 越 小 表明 干旱 程度 越 大 ,本 人 研究 采用 
Kogen 在 文献 * 中 报道 的 0.35 作 为 闷 值 判断 干旱 是 
否 发 生 , 即 VCI< 0.35 时 干旱 发 生 ,VCI>0.35 时 干旱 
不 发 生 。 基 于 干旱 发 生 的 栅 格 数 , 计 算 区 域 干旱 面 
积 占 比 。 

2.2.2 和 干旱 风险 评估 方法 干旱 灾害 风险 由 致 灾 因 
TIEKE (Pay) 、 承 载体 脆弱 性 (Es) 孕 灾 环境 暴露 
度 (有 is) 人 类 防 灾 减 灾 能 力 (4u)4 部 分 构成 ”2 ,本 
研究 假设 各 种 因素 对 干旱 风险 的 贡献 相等 ,干旱 灾 
害 风 险 指 数 R, 由 下 式 计算 : 

R,, = Stand(P,,) X E,, XV,, x( 2 A,,) (2) 

本 研究 基于 干旱 发 生 频率 进行 干旱 危险 性 评 
佑 ,计算 方法 宛如 下 : 


met (3) 


式 中 :Ne 为 总 年 数 ;me 为 发 生 干旱 的 年 数 ;Pu 高 表 
明 干 旱 危 险 性 高 , 低 则 表明 干旱 危险 性 低 。 

考虑 到 在 巴基斯坦 所 有 作物 中 ,小麦 种 植 范围 
最 广 ,而 小 麦 产量 与 2 月 NDVI 密 切 相 关 ” ,本 研 
究 利 用 2 月 NDVI 获取 Ps; 并 利用 如 下 标准 化 函数 
Stand 对 Ps 等 变量 进行 标准 化 。 


Stand(X)=0.5+0.5x 


X =A un 

X, Xn. v 

RP Ls X OU BERE DUREE X BY MURUS TÉ. 
FBLA WEBS PEU, PEU SAT : 


E,, =Stand(DA,,,) x Stand(Goil,,,) (5) 
š 0 = Omin 
oily, E B x ES O min d 


UP DAS ACHE 8 LE Bi (96) ; Soil, 为 标准 化 的 土 
壤 田 间 持 水 量 ;9 为 土壤 田间 持 水 量 (%) ,根据 土壤 


质地 资料 获取 ,0 和 0 分 别 为 研究 区 最 大 和 最 
小 田间 持 水 量 (%)。 
孕 灾 环境 暴露 度 (V,) 计 算 方法 如 下 : 
Viy = Stand(R nwa) X Stand(PD ,,,) x Stand(Yield,,) (7) 
PD - PD,, 
PD,, = PD. PD. PD (8) 
Yield - Yield, 
ed. Yield. (9) 
式 中 :Rie 为 作物 种 植 百 分 比 (%) ;PD 为 标准 化 
的 人 口 密度 (人 :km”);PD 为 人 口 密度 (人 :km”); 
PD 、PD% 分 别 为 研究 区 最 小 、 最 大 人 口 密度 (人 … 
km?) ; Yield, 为 标准 化 的 综合 粮食 产量 (kg: hm?) ; 
Yield 为 某 地 区 综合 粮食 产量 (kg* hm?) ; Yield, 
Yieldw 分 别 为 研究 区 最 小 、 最 大 综合 粮食 产量 (kg- 
hm?), 
人 类 防 灾 减 灾 能 力 4w 计 算 方法 如 下 : 
A,, =Stand(R,,) (10) 
式 中 :R, 代 表 灌 溉 面积 比例 (% ) o 


3 结果 与 分 析 


3.1 巴基斯坦 月 尺度 干旱 范围 及 影响 因子 

如 图 4 所 示 , 巴 基 斯 坦 在 1982 一 2013 年 ,1 月 和 
9 一 12 月 干旱 面积 比例 大 于 40%,12 月 干 星 面 积 比 
例 为 全 年 最 高 (面积 比例 45%);5 一 8 月 干旱 面积 比 
例 较 小 ,6 月 干旱 面积 全 年 最 小 ,比例 小 于 30%。 通 
过 偏 相关 系数 分 析 , 各 月 干旱 面积 比例 与 各 月 气 
温 、 降 水 的 偏 相 关系 数 分别 为 0.69 .0.74, 分 别 通过 
了 5%、1% 水 平 显著 性 检验 ,说 明 干 旱 范 围 同 时 受气 
温和 降水 的 影响 ,但 降水 的 影响 程度 大 于 气温 。 

1982 一 2013 年 ,巴基斯坦 各 月 干旱 面积 比例 变 


首届 


2 3 4 5 6 7 8 9 10 11 12 
月 份 


Yield, 


多 年 平均 干旱 面积 比例 /% 


图 4 巴基斯坦 1982 一 2013 多 年 平均 各 月 干旱 面积 比例 
Fig. 4 Average proportion of dry areas in each month in 
Pakistan during 1982—2013 
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化 不 一 (图 5)。5 一 6 月 干旱 面积 比例 呈 1% 水 平 显 
车 上 升 ,5 月 上 升 速率 最 大 ,6 月 次 之 。 造 成 该 现象 
的 主要 原因 是 温度 显著 上 升 ,但 降水 呈 下 降 趋势 ， 
温度 升 高 和 降水 减少 导致 5 一 6 月 干旱 影响 范围 扩 


F 旱 特征 及 其 风险 评估 


32 巴基斯坦 干旱 频率 
由 图 7 和 土地 覆盖 空间 分 布 图 可 以 看 出 ， 
1982—2013 年 巴基斯坦 西部 山区 林地 .草地 9 月 干 


大 。 其 中 6 月 是 干旱 面积 比例 最 小 月 (图 4) ,干旱 面 
积 的 扩大 意味 着 6 月 干旱 状况 加 重 ,这 将 导致 巴 基 
斯 坦 旱 灾 加 剧 。1 月 和 2 月 干旱 面积 分 别 在 1% 和 
5% 水 平 上 呈 显 著 下 降 , 其 中 1 月 下 降 速 率 较 大 ;其 
余 月 干旱 面积 比例 变化 不 显著 。 

为 分 析 1 月 干旱 面积 比例 变化 原因 ,分 析 了 1 
月 VCI 的 变化 。 从 图 6 可 以 看 出 ,1982 一 2013 年 巴 
基 斯 坦 1 月 VCI 呈 明显 上 升 趋势 的 地 区 主要 位 于 
Punjab 省 ,说 明 Punjab 省 干旱 程度 减 小 , 且 干 旱 程度 
减 小 的 面积 呈 扩 大 趋势 。 因 此 ,巴基斯坦 1 月 干旱 
面积 显著 减少 主要 源 于 Punjab 省 干旱 面积 的 缩 
减 。 根 据 已 有 报道 ”下 ,1980 年 以 来 ,Punjab 省 冬季 
地 表 供 水 量 呈 减少 趋势 ,灌溉 所 用 地 下 水 开采 量 呈 
上 升 趋势 ” ,考虑 到 1 月 之 前 有 2 次 灌溉 ,说 明 地 下 
水 开采 量 的 增加 使 得 Punjab 省 1 月 干旱 面积 显著 


旱 频 率 较 高 ,将 影响 当地 生态 环境 ;1 月 和 9 一 12 月 ， 
东部 平原 降水 较 少 ,干旱 频率 较 高 , 且 这 些 月 多 年 
平均 干旱 面积 比例 较 大 ,考虑 到 东部 平原 作物 种 植 
面积 比例 大 ,干旱 将 影响 该 地 区 作物 产量 。 此 外 ,9 
月 干旱 频率 较 高 的 范围 最 广 ,主要 分 布 于 西部 山区 
向 东部 平原 区 过 渡 带 ,将 导致 东部 平原 区 西 侧 生 态 
环境 脆弱 ,影响 东部 作物 生长 。5 一 8 月 为 巴基斯坦 
的 夏季 , 受 季 风 降 水 影响 ,全 境 干旱 频率 较 低 , 有 利 
于 当地 社会 经 济 生产 活动 。 
3.3 巴基斯坦 干旱 风险 评估 

从 图 8 可 以 看 出 ,1982 一 2013 年 巴基斯坦 干旱 
风险 较 高 的 地 区 位 于 北部 山区 向 平原 过 渡 的 地 带 
以 及 东南 部 自然 植被 覆盖 地 区 ;干旱 风险 较 低 的 地 
区 主要 位 于 地 形 平坦 的 中 东部 河谷 地 区 以 及 北部 
山区 。 对 整个 巴基斯坦 而 言 ,干旱 风险 空间 格局 与 
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图 5 1982 一 2013 年 各 月 干旱 面积 比例 变化 
Fig. 5 Variation of proportion of dry areas in Pakistan in each month during 1982—2013 
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图 6 巴基斯坦 1 月 多 年 平均 植被 状态 指数 (VCI) 


Fig.6 Pakistan's multi-year average vegetation condition index (VCI) in January 
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图 7 基于 VCI 的 巴基斯坦 1982 一 2013 年 各 月 干旱 频率 
Fig. 7 Drought frequency of each month in Pakistan during 1982—2013 based on VCI 


干旱 脆弱 性 、 危 险 性 空间 格局 相似 ,特别 是 干旱 脆 。 ”西部 的 Balochistan 省 ,干旱 风险 较 低 , 但 干旱 暴露 度 
弱 性 。 但 干旱 风险 空间 格局 与 干旱 暴露 度 空间 格 。” 高 。 根 据 偏 最 小 二 乘法 获得 的 各 因素 重要 性 指标 
局 相差 较 大 ,特别 是 东部 的 Punjab 和 Sindh 省 ,以 及 (Variable importance index, VIP) ,巴基斯坦 干旱 风 
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二 ”标准 化 旱灾 暴露 度 
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F 旱 特征 及 其 风险 评估 
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图 8 1982—2013 年 巴基斯坦 干旱 风险 
Fig. 8 Drought risk in Pakistan during 1982—2013 


险 空间 格局 主要 由 干旱 频率 、 作 物产 量 、 大 牲畜 比 
例 、 土 壤 持 水 能 力 (VIP>1) 决 定 ( 图 9) ,其 中 干旱 频 
率 对 干旱 风险 的 影响 最 大 。 


干旱 脆弱 性 均 较 高 (图 8a~b), 需 适度 放牧 或 减少 大 
牲畜 数量 。 在 东南 部 自然 植被 区 ,土壤 质地 为 砂 土 
SUR, E ,土壤 持 水 能 力 较 低 , 导 致 干旱 脆弱 性 较 


分 析 表 明 ,在 北部 山区 ,自然 植被 覆盖 度 较 高 ， 
作物 产量 低 , 导 致 干旱 风险 性 和 干旱 暴露 度 均 较 低 
(图 8a\c); 但 该 地 区 土壤 质地 为 壤 砂 ,导致 干旱 脆 
弱 性 较 高 , 需 开 展 水 土 保持 工作 。 在 北部 山区 向 平 
原 过 渡 的 地 区 ,大 牲畜 比例 较 高 ,导致 干旱 风险 和 


VIP 
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ORO oi 


全 & > K & 
SD Que Aw uo Sa 
xe Fl 9 这 ge 

x 


干旱 风险 要 素 


注 :VIP 为 各 要 素 重 要 性 指标 。 

图 9 基于 偏 最 小 二 乘法 的 巴基斯坦 干旱 风险 要 素 
重要 性 指标 (VIP) 

Fig.9 Variable importance index (VIP) of drought risk 


elements in Pakistan based on partial least square method 


高 ,同时 该 地 区 干旱 频率 较 高 ,导致 干旱 危险 性 较 
高 , 较 高 的 脆弱 性 和 和 危险 性 使 该 地 区 干旱 风险 较 
高 ,为 降低 该 地 区 干旱 风险 并 改善 生态 环境 质量 ， 
需 增 加 生态 环境 配水 ,开展 水 土 保持 工作 。 在 各 因 
素 的 综合 作用 下 ,中 西部 干旱 风险 处 于 中 间 水 平 , 考 
虑 到 该 地 区 以 自然 植被 为 主 ,土壤 质地 为 砂 士 或 壤 
砂 土 ,土壤 持 水 能 力 较 低 , 也 需 进行 水 土 保持 ”””。 
在 地 形 平坦 的 中 东部 河谷 地 区 ,如 Punjab 省 中 东 
部 .Sindh 省 中 部 ,作物 产量 较 高 ,导致 干旱 暴露 度 较 
高 ;风险 和 干旱 危险 性 都 较 低 (图 8a、d) ,因为 该 地 
区 靠 开 采 地 下 水 维持 较 高 的 灌溉 面积 比例 ,但 是 开 
采 地 下 水 将 消耗 大 量 的 能 量 , 超 采 的 地 下 水 将 引起 
盐 渍 化 ,导致 环境 退化 。 印 度 河 是 巴基斯坦 重要 的 
地 表 灌 溉 水 源 ,在 气候 变化 和 人 类 活动 影响 下 ,位 
于 印度 河 下 游 的 巴基斯坦 中 东部 河谷 地 区 获得 的 
水 资源 量 逐 渐 减 少 , 从 可 持续 发 展 的 角度 ,相关 部 
门 需 加 强 印度 河流 域 尺度 的 水 资源 管理 ,改善 河道 
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沿岸 各 用 水 机 构 对 地 表 水 的 管理 ,使 目前 分 配 不 到 
地 表 水 用 户 的 农业 生产 力 上 升 ,从 而 提高 水 土 .能 
量 的 利用 效率 ,以 应 对 下 游 地 下 水 日 益 减 少将 引发 
的 环境 退化 及 高 干旱 风险 ”。 


4 讨论 


当前 大 部 分 研究 只 是 针对 巴基斯坦 干旱 较为 
严重 的 部 分 区 域 展开 干旱 特征 分 析 , 没 有 对 巴 基 斯 
坦 全 国 展开 干旱 分 析 , 且 没有 对 巴基斯坦 展开 干旱 
灾害 风险 分 析 。 如 Shah 等 中 和 Adnan 等 '" 均 对 巴 基 
斯 坦 Sindh 省 干旱 状况 展开 研究 ,研究 结果 均 表 明 
Sindh 省 干旱 频率 高 ,干旱 状况 严重 ,尤其 是 Sindh 省 
东南 部 干旱 程度 较 高 ,这 与 本 文 的 干旱 风险 空间 格 
局 是 一 致 的 。 但 由 于 数据 方法 和 研究 区 范围 的 差 
别 ,研究 结果 也 存在 一 定 的 差异 。Shah 等 ”认为 由 
F Sindh 省 高 度 依赖 季风 降水 ,季风 降水 的 不 足 导 
致 干 旱 的 发 生 , 且 在 季风 季节 ,北部 降水 比 南部 少 ， 
北部 比 南部 更 容易 遭受 干旱 ,因此 Sindh 省 北部 干 
旱 频 率 高 于 南部 。 但 从 本 文 各 月 干旱 频率 空间 分 
布 ( 图 7) 和 王 旱 风险 分 布 ( 图 8) 来 看 ,Sindh 省 南北 
于 旱 频 率 差别 不 大 ,这 可 能 与 本 文 充分 考虑 到 灌溉 
影响 有 关 。 值 得 一 提 的 是 ,Adnan 等 "在 研究 气温 、 
降水 和 干 星 的 相关 性 时 提 到 了 土壤 持 水 性 和 地 表 
水 灌溉 对 巴基斯坦 农业 的 影响 ,但 并 没有 展开 分 
析 。 本 文采 用 时 间 和 空间 尺度 较 大 的 CIMMS NDVI 
长 时 间 序 列 数据 集 , 将 遥感 数据 与 社会 经 济 资料 相 
结合 ,综合 考虑 干旱 频率 作物 产量 、 大 牲畜 比例 、 
土壤 持 水 能 力 灌溉 等 因素 对 干旱 的 影响 ,对 巴 基 
斯 坦 全 国 1982 一 2013 年 展开 长 时 间 尺 度 的 干旱 特 
征 及 风险 研究 ,分析 了 巴基斯坦 全 国 和 干旱 及 其 风险 
影响 因素 ,为 巴基斯坦 干旱 监测 以 及 灾害 管理 机 构 
提供 了 有 用 的 参考 。 

本 研究 收集 的 统计 资料 有 限 , 且 由 于 实测 气象 
数据 缺乏 ,所 用 再 分 析 资 料 以 及 干旱 指数 在 研究 中 
存在 的 问题 ,将 在 一 定 程度 上 影响 研究 结果 的 准确 
性 。 因 此 该 研究 仅 在 现 有 数据 的 基础 上 ,对 巴 基 斯 
坦 长 时 间 尺 度 的 干旱 风险 进行 评估 ,为 巴基斯坦 自 
然 灾害 风险 防 控 和 减灾 对 策 提供 依据 。 在 未 来 的 
研究 中 ,考虑 利用 时 间 序 列 更 加 完善 的 统计 实测 
数据 等 ,研究 巴基斯坦 干旱 风险 时 空 变化 ,为 未 来 
气候 变化 下 巴基斯坦 干旱 风险 管理 提供 更 加 完备 
的 科技 支撑 。 


5 结论 


< 


为 了 加 强 巴 基 斯 坦 干 旱灾 害 风 险 管理 ,促进 当 
地 社会 经 济 发 展 ,本 文 利 用 1982 一 2013 年 GIMMS 
NDVI 数 据 集 ,计算 植被 状态 指数 (VCI) ,采用 VCI 研 
究 区 域 干旱 时 空 变化 特征 ,并 通过 VCI 获 取 干 量 危 
险 性 ,利用 作物 产量 、 作 物种 植 面积 等 统计 资料 获 
KTERE .脆弱 性 等 指标 ,进而 基于 自然 灾害 
风险 理论 ,对 巴基斯坦 干旱 灾害 风险 进行 空间 分 
析 , 得 到 以 下 结论 : 

(1) 1982 一 2013 年 ,受气 温和 降水 的 影响 , 巴 基 
斯 坦 大 部 分 月 份 干旱 面积 比例 均 在 30% 以 上 ,干旱 
范围 大 。 各 月 干旱 面积 比例 变化 不 一 ,5 一 6 月 干旱 
面积 比例 显著 上 升 ,1 一 2 月 干旱 面积 比例 显著 下 
降 ,其 余 月 份 变化 不 显著 。 由 于 开采 地 下 水 进行 灌 
溉 ,Punjab 省 1 月 于 量程 度 减弱 ,成 为 该 月 巴基斯坦 
干旱 面积 比例 下 降 的 重要 原因 。 

(2) 根据 干旱 频率 表征 的 干旱 危险 性 ,西部 山 
区 林地 、 草 地 9 月 干旱 危险 性 较 高 ,东部 平原 1 月 、 
9 一 12 月 干旱 危险 性 较 高 ,这 些 月 同时 为 多 年 平均 
干旱 面积 比例 较 大 的 月 。 

(3) 由 危险 性 .脆弱 性 .暴露 度 和 人 类 防 灾 减 灾 

得 到 干旱 风险 评估 结果 表明 ,巴基斯坦 干旱 风 
仿 较 高 的 地 区 位 于 北部 山区 向 平原 的 过 渡 地 带 以 
及 东南 部 自然 植被 覆盖 的 地 区 ,干旱 风险 间 格 局 主 
要 受 干旱 频率 、 作 物产 量 、 大 牲畜 比 例 .土壤 持 水 能 
力 影响 。 

GIMMS NDVI3g 是 已 知 最 长 时 间 序 列 的 NDVI 
数据 集 * 外 , 且 具 有 稳定 性 ,对 干旱 监测 的 连续 性 较 
好 。 本 文 将 该 遥感 资料 与 社会 经 济 资料 相 结合 ,从 
灾害 风险 的 角度 分 析 和 评估 了 巴基斯坦 的 干旱 灾 
害 ,较为 全 面 的 反映 了 自然 和 人 为 因素 作用 下 的 巴 
基 斯 坦 干旱 风险 分 布 。 本 研究 是 目前 针对 巴 基 斯 
坦 全 境 最 长 时 间 跨 度 的 干旱 特征 研究 , 且 评 佑 了 当 
前 学 者 尚未 涉及 的 巴基斯坦 干旱 风险 ,为 巴基斯坦 
干旱 灾害 风险 管理 提供 了 有 益 参 考 。 
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Drought characteristics and its risk assessment across Pakistan 
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Abstract: Pakistan is in South Asia (60?—80?N, 23°-37°E) with a total area of 88x10* km’. The southeastern part 
of Pakistan is mountainous and hilly with a complex topography and different climate types. In summer, this area 
is controlled by subtropical high pressure, which results in hot and dry weather with a small amount of 
precipitation. The average annual rainfall is less than 250 mm, and a quarter of the area receives less than 120 
mm. Because of the high summer temperatures, Pakistan has been greatly affected by long-term drought, which 
has severely restricted its social and economic development. In this study, we used the Global Inventory 
Modeling and Mapping Studies Normalized Difference Vegetation Index 3" Generation (GIMMS NDVI3g) 
dataset from 1982 to 2013 to construct the vegetation condition index (VCI). We then used the VCI to analyze the 
range and frequency of drought in Pakistan for the past 32 years. The drought frequency was used to characterize 
the risk of drought based on sources from the Pakistan Statistics Bureau and reports in the literature. Information 
on livestock, soil, crops, and irrigation were used to evaluate the drought vulnerability, exposure, and disaster 
prevention and mitigation capabilities of Pakistan. Additionally, ERA-Interim data that had been verified by 10 
meteorological stations in Pakistan were used for temperature and precipitation reanalysis. The influencing 
factors of the drought range were determined using partial correlation coefficients, and the partial least squares 
method was used to identify the dominant influencing factors for the spatial patterns of regional drought risk. The 
following results were obtained. (1) During the study period, the multiyear average drought range differed for 
each month. The drought range was larger from January to September than from September to December (area 
ratio was greater than 40%). The drought range was smaller from May to August (area ratio was less than 30%). 
The increase in temperature and decrease in precipitation from May to June caused the arid area to increase. (2) 
In the study area, the drought frequency was low from May to August. The frequency was higher in September 
for the woodlands and grasslands of the western mountainous areas. The frequency was higher in January and 
September than in December for the eastern plains. (3) The spatial pattern of drought risk was mainly determined 
by the drought frequency, crop yield, proportion of large livestock, and water storage capacity of the soil. Among 
these, the drought frequency had the highest impact on the regional drought risk. The livestock ratio and the water 
storage capacity of the soil indicated a higher drought risk in the transition zone from mountains to plains in the 
northern area and the southeastern area covered by natural vegetation. Irrigation reduced the drought risk in the 
valley area with flat terrain. Additionally, the natural vegetation in the northern mountainous area reduced the 
drought risk as well. In this study, we analyzed the long-term drought characteristics in Pakistan and evaluated 
the drought risk due to natural and human factors. This research can provide useful information for drought and 
disaster risk management in Pakistan. 


Key words: drought; risk assessment; normalized difference vegetation index(NDVI); vegetation condition in- 
dex(VCI); Pakistan 


